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Using 22.7 million BB events recorded with the BABAR detector, the inclusive branching fractions 

for the production of J/ifi, ip(2S) and Xc in B decays are presented. Combining the charmonium 

43 ■ state with & K± , K° , K*± , K*° or 7r , B decays are reconstructed exlusively and branching fractions 

are determined. A preliminary study is also presented for the B — > r/ c K decay mode. 
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1 Introduction 



Reconstruction and study of charmonium mesons in B decays is a crucial component of the mea- 
surement of time-dependent CP-violating asymmetries [lj. 

The analyses described in the following paper are based on a sample of 20.7 fb _1 collected with 
BABAR at the Y(4S) resonance with an additional 2.6 fb _1 collected below the BB threshold. A 
determination of the B meson branching fractions depends upon an accurate measurement of the 
number of B mesons in the data sample. The number of BB events is determined by comparing the 
rate of multi-hadron events in data collected both on and off resonance. The continuum contribution 
to the on-resonance sample is estimated by rescaling the number of off-resonance hadronic events 
by the ratio of the number of observed fJ- + fJ>~ events in the two samples. This procedure yields a 
total of 22.72 ± 0.36 million BB events. 



2 Inclusive decays of B to states containing Charmonium 

J/ip candidates are selected by requiring two identified leptons of opposite charge. Electrons are 
identified based on the ratio E/p of the energy deposited in the calorimeter to the measured 
momentum from tracking information, the shape of the calorimetric cluster and the ionization in the 
tracking detectors. Muons are identified by requiring a minimum ionizing signal in the calorimeter. 
In addition, the shape and penetration of the distribution of hits in the instrumented flux return 
are used. The number of J/ip events is determined by fitting the invariant mass distribution to 
a probability density function obtained from a simulation including contribution from both final 
state radiation and bremsstrahlung. The fit yields 15739 ± 177 J/ip — > e + e~ and 13683 ± 154 J/ip 
— > signal events (Figure ||). 
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Figure 1: Invariant mass distributions for inclusive J/ip production in B decays 

The ip(2S) candidates are reconstructed in the tp(2S) — > l + l~ and ip(2S) — > J/ipir + -K~ decays. 
In the latter case the signal yield is extracted by a fit to the mass difference between the ip(2S) and 
the J/ip reconstructed candidates. We find 552 ± 50 decays to e~e~ , 437 ± 44 decays to 
405 ± 37 decays to J/ip(e + e~)ir + ir~ and 400 ± 34 decays to J/V'()U + A t_ ) 7r+7r_ - 

The Xci and Xc2 candidates are selected by combining a reconstructed J/ip with a photon. The 



signal yield is determined by fitting the mass difference between the \ c and the J/ifj candidates. 
The shape of the signal is extracted from Monte Carlo, with the mass difference between the Xci and 
Xc2 peaks fixed to the PDG value Q. The fit gives 476 ± 71 Xci and 86 db 59 Xc2 candidates for the 
J/ip — > e + e~ decay and 545 ± 60 Xci and 104 ± 56 Xc2 candidates for the J/ip — > decay. 



Mode 


Br(xlO^) 


B — > J/tp X 
B ip(2s) X 
B - XdX 

B - Xc2^ 


1.044 ± 0.013 ± 0.028 
0.275 ±0.020 ±0.029 
0.378 ± 0.034 ± 0.026 
< 0.21 (90% CL) 



Table 1: Measured branching fractions for inclusive charmonium production in B decays. 

Values for the branching fractions are extracted from the yields ||. A 90% CL limit is set on 
the B decay to Xc2 ■ The results for the inclusive production of J/ip, ip(2S) and Xc i n B decays are 
summarized in Table |], where the first uncertainty is statistical error and the second is systematics. 

3 Exclusive decays of B to Charmonium 

The reconstruction of exclusive decay modes containing charmonium presents a very low back- 
ground in most of the channels. For this reason the lepton identification criteria are loosened for 
one of the two J/ijj decay products. As in the inclusive analysis, ip{2S) candidates are reconstructed 
by their decay to e~e~, and J/ip ir + i:~ . The Xc candidates are selected through their decay 

to J/?/>7- 

The charmonium states are selected in a window around their expected mass|Q] for the decays 
to leptons. In the decays to J/ip states the mass difference distribution between the charmonium 
candidate and the reconstructed J/ip is used instead. 

Selected candidates are then paired with a K + , K® (either ir + ir~ or 7r 7r° ), K* + (either K + ir° 
or 7r + ), K*° (either K + ir~ or K® ir ), 7r° or Kl to form a B candidate. The two most significant 
observables used to identify the signal are the difference AE in the center-of-mass frame between the 
reconstructed B energy and half the nominally available energy, \/s/2, and the energy-substituted 

mass, rriEs = ys /4 — P B 2 , where Pg is the center-of-mass momentum of the B candidate. A 
sample of these distributions is given in Figure [2] for J/ipKs events. 

In the case of multiple candidates per event, only the candidate with the smallest \AE\ is 
selected. For all modes except B — > J/ipKi and B — > J/ipK*, the number of signal events 
is determined from the observed number of events in the (AE,ities) region after background 
subtraction. In addition to the usual combinatorial component, the background distribution shows 
an excess of events in the signal region. The combinatorial contribution is estimated by using 
an ARGUS function in the fit to the ities distribution. The peaking background component is 
obtained from simulation of inclusive B decays to charmonium, after removing the signal events. 

The signal yields for the K*° and K* + modes are determined simultaneously from a likelihood 
fit, which is needed to account for the cross-feed between the K* decay channels. 

A different technique is used for the B — > J/i/jKl decay mode. In this case only the Kl direction 
is measured with information from the calorimeter and the instrumented flux return. Given this 
direction and the reconstructed charmonium candidate, the Kl energy is extracted by using the B 
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Figure 2: Example of the AS (left) and mgg (right) distributions for the decay B° — ► J/ipK® 

mass as a constraint. To eliminate cross- feed from other decay modes, a veto has been introduced 
for events which have been selected already in other exclusive modes. This procedure yields a 
purity of about 50%. Due to this method, the thes distribution cannot be used to determine the 
signal yield. The AEj^ K o distribution is used in a log-likelihood fit. The shapes of the signal 
and inclusive charmonium background components are taken from Monte Carlo simulations. The 
shape of the non-charmonium background component is taken from an ARGUS fit to the AE K o 
distribution for events in the J/ip mass sideband. After the background subtraction, this channel 
gives a signal yield of 183 ± 14 events (Figure ||, left plot). 




Figure 3: (Left plot) AE distribution for B — ► J/i^Kl decays. Points are data, solid line is the 
Monte Carlo simulation. The three components are respectively signal, background events which 
include a real J/tp, and non-J/^ background. (Right plot) preliminary study of B decay to rj c K. 
In the above plot the r] c candidate is reconstructed through KsKir decays. 



The information on the signal yields and the measured branching fractions for all presented 



Mode 


Yield 


Br(xl0^ 4 ) 


B u 


— ^ J/ipK v 


K" s -> 7T + 7T- 


265.5 ± 2.9 


8.5 ±0.5 ±0.6 






#g -> 7T 7r° 


62.5 ±3.8 


9.6 ±1.5 ±0.7 








183 ± 14 


6.8 ±0.8 ±0.8 






All 




8.3 ±0.4 ±0.5 


B+ 


-> J/iPK+ 




1109 ±4 


10.1 ±0.3 ±0.5 


B° 






13.6 ±0.9 


0.20 ±0.06 ±0.02 


B" 






594 ± 8.5 


12.4 ±0.5 ±0.9 


B+ 


-> J/^K*+ 




377.4 ± 16.9 


13.7 ± 0.9 ± 1.1 


B° 


-» V(2S)iT 




56.0 ±3.4 


6.8 ±1.0 ±1.1 


B+ 


-> ip(2S)K + 




207.3 ±6.2 


6.3 ±0.5 ±0.8 


B 1 ' 


- Xci^° 




26.1 ± 2.5 


5.4 ± 1.4 ± 1.1 


B+ 


- Xd^ + 




145.1 ± 7.2 


7.5 ±0.8 ±0.8 


B° 






32.6 ±6.0 


4.8 ±1.4 ±0.9 


B° 


— » J/if}ir + ir 




29.1 ±9.4 


0.46 ± 0.11 ± 0.08 



Table 2: Summary of the signal yields and extracted branching fractions for the different B decays 
reconstructed into exclusive charmonium final states by the BABAR analyses described in this paper. 

exclusive modes [|j is summarized in Table [2], while Figure || shows a comparison of the new pre- 
liminary BABAR results to the current PDG values. 
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Figure 4: Comparison of the preliminary results from the BABAR summarized in the previous table 
and the published PDG values. 



In addition to the analyses described above, a very preliminary study has been recently per- 
formed by BABAR, on exclusive reconstruction of B decays to r/ c K modes. The r] c has been studied 
in the KsK + ir~ , K + K~-k° and K + K~ K + K~ {(fxj) plus non-resonant) decay modes. The energy- 
substituted mass distribution for the decay i] c — > KsK + ir~ is presented as an example in Figure |||, 
right plot. No value has been extracted yet for the branching fractions in these decay modes. 

4 Summary 

Using 22.7 million BB events recorded by the BABAR detector, the inclusive branching fractions for 
the production of J/tp , ip{^S) and Xc are presented. Combining the charmonium state with either 
a , K°, K*^ , K*° or 7T°, B decays are reconstructed exclusively and their branching fractions 
are determined. 
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